SMC-TR-99-12 


AEROSPACE  REPORT  NO. 
TR-98(8555)-2 


Semiclassical  Random  Electrodynamics: 
Spontaneous  Emission  and  the  Lamb  Shift 


15  April  1999 


Prepared  by 

J.  C.  CAMPARO 
Electronics  Technology  Center 
Technology  Operations 


Prepared  for 

SPACE  AND  MISSILE  SYSTEMS  CENTER 
AIR  FORCE  MATERIEL  COMMAND 
2430  E.  El  Segundo  Boulevard 
Los  Angeles  Air  Force  Base,  CA  90245 


19990521  1/2 


Engineering  and  Technology  Group 


THE  AEROSPACE 

CORPORATION 

egundo,  California 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED 


QUALITY 


Ws*bctbd 


4 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering 
and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite 
1204,  Arlington,  VA  22202*4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704*0188),  Washington,  DC  20503. 

1  AGENCY  USE  ONLY  {Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

15  April  1999 

4.  TITLE  AND  SUBTITLE 

Semiclassical  Random  Electrodynamics: 

Spontaneous  Emission  and  the  Lamb  Shift 

5.  FUNDING  NUMBERS 

F04701 -93-C-0094 

6.  AUTHOR(S) 

J.  C.  camparo 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

The  Aerospace  Corporation 

Technology  Operations 

El  Segundo,  CA  90245-4691 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

TR-98(8555)-2 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Space  and  Missile  Systems  Center 

Air  Force  Materiel  Command 

2430  E.  El  Segundo  Boulevard 

Los  Angeles  Air  Force  Base,  CA  90245 

io.  spons6ring/monitoring 

AGENCY  REPORT  NUMBER 

SMC-TR-99- 1 2 

11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

It  is  often  remarked  that  an  explanation  of  spontaneous  emission  and  the  Lamb  shift  requires  quan¬ 
tization  of  the  electromagnetic  field.  Here,  these  two  quantities  are  derived  in  a  semiclassical  form¬ 
alism  by  use  of  second-order  perturbation  theory.  The  purpose  of  this  report  is  not  to  argue  the  val¬ 
idity  of  QED  but  rather  to  develop  a  semiclassical  approximation  to  QED  that  may  nonetheless  have 
certain  computational  advantages  over  QED.  To  this  end,  the  vacuum  of  QED  is  simulated  with  a 
classical  zero-point  field  (ZPF),  and  as  a  consequence,  the  resulting  theory  is  entitled  semiclassical 
random  electrodynamics  (SRED).  In  the  theory,  the  atom  is  coupled  to  the  ZPF  and  to  its  own 
radiation-reaction  field  through  an  electric  dipole  interaction.  These  two  interactions  add  to  produce 
exponential  decay  of  excited  states  while  they  cancel  each  other  to  prevent  spontaneous  excitation 
of  the  ground  state;  the  Lamb  shift  appears  in  the  theory  as  an  ac  Stark  shift  induced  by  the  ZPF. 

The  spontaneous  decay  rate  of  an  excited-state  derived  in  SRED  is  equal  to  the  Einstein  A  coeffi¬ 
cient  for  that  state,  and  the  Lamb  shift  agrees  with  that  of  nonrelativistic  QED.  Moreover,  SRED  is 
shown  to  be  useful  for  the  numerical  simulation  of  spontaneous  decay. 

14.  SUBJECT  TERMS 

Lamb  shift,  quantum  electrodynamics,  semiclassical  random  electrodynamics 

15.  NUMBER  OF  PAGES 

9 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED  UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

NSN  7540-01  -280-5500  Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  Z39-1 8 


298-102 


Contents 


1.  INTRODUCTION .  1 

2.  RADIATION  REACTION .  1 

3.  CLASSICAL  ZERO-POINT  FIELD .  2 

4.  PERTURBATION  THEORY .  3 

A.  Evolution  of  the  Wave  Function .  3 

B.  Wave-Function  Evolution  Owing  to  Radiation  Reaction .  4 

C.  Wave-Function  Evolution  Owing  to  the  Zero-Point  Field .  4 

D.  Spontaneous  Decay .  6 

E.  Lamb  Shift .  6 

5.  NUMERICAL  SIMULATION  OF  SPONTANEOUS  DECAY .  6 

6.  SUMMARY .  8 

REFERENCES  AND  NOTES .  8 

Figures 

2 

1 .  Numerical  simulation  of  the  sodium  3  P,/2  state’s  spontsaneous  decay  using  SRED .  8 


iii 


1.  INTRODUCTION 

Though  few  would  argue  that  quantum  electrodynamics 
is  an  accurate  description  of  nature,  the  semiclassical 
theory  of  field-atom  interactions  still  holds  an  important 
place  in  modern  physics.  One  need  only  consider  the 
Bloch  equations  and  their  utility  in  quantum  optics1,2  and 
magnetic  resonance3  to  appreciate  the  vitality  of  semi¬ 
classical  electrodynamics.  Notwithstanding  this  utility, 
a  failing  of  traditional  semiclassical  electrodynamics4  has 
been  its  inability  to  predict  accurately  spontaneous  emis¬ 
sion  from  first  principles.5,6  There  is  nonetheless  a  feel¬ 
ing  among  many  that  semiclassical  electrodynamics 
ought  to  be  more  successful  in  this  regard;7  after  all,  QED 
indicates  that  vacuum  fluctuations  and  radiation  reaction 
are  inextricably  linked  (and  therefore  that  spontaneous 
decay  can  be  interpreted  as  an  effect  driven  by  radiation 
reaction8)  and  radiation  reaction  is  a  feature  of  classical 
electrodynamics.9 

Beyond  the  issue  of  classical  field  theory’s  intrinsic  ca¬ 
pabilities,  there  are  pragmatic  reasons  for  pursuing  a 
semiclassical  description  of  spontaneous  emission.  Typi¬ 
cally,  spontaneous  decay  is  incorporated  into  the  semi¬ 
classical  theories  of  quantum  optics  through  an  exponen¬ 
tial  decay  rate  for  atomic  density  matrix  elements. 
However,  there  are  processes  such  as  quantum  jumps10 
and  radiative  decay  in  photonic  bandgaps11  for  which  this 
approach  is  inappropriate.  Although  valid  and  success¬ 
ful  theoretical  methods  exist  for  describing  these  pro¬ 
cesses,  their  semiclassical  interpretation  might  neverthe¬ 
less  have  utility.  In  this  regard,  the  present  work  is 
meant  as  a  step  toward  extending  the  reach  of  semiclas¬ 
sical  methods  to  a  broader  class  of  radiative  phenomena. 

To  this  end,  it  is  worth  noting  that  no  matter  how  one 
interprets  the  cause  of  spontaneous  decay  in  QED,  the 
vacuum  plays  a  significant  role.12,13  Consequently,  if  a 
semiclassical  theory  is  to  be  successful  in  describing  this 
phenomenon,  then  vacuum  fluctuations  should  somehow 
appear  in  the  theory.14  Here  the  vacuum  of  QED  is 
simulated  with  a  classical  zero-point  field  (ZPF),15  and 
the  consequences  of  this  ZPF  in  semiclassical  electrody¬ 
namics  are  explored.  For  clarity,  the  resulting  theory  is 
titled  semiclassical  random  electrodynamics  (SRED).  I 
wish  to  state  explicitly  that  I  make  no  claims  for  SRED  as 
a  rival  to  QED;  the  present  experimental  situation  pro¬ 
vides  no  imperative  for  such  a  rival.  If  anything,  low-to- 
high-energy  experiments  argue  forcefully  for  the  validity 
of  quantized  fields.  Rather,  SRED  should  be  viewed  as  a 
semiclassical  approximation  to  QED,  which  may  never¬ 
theless  have  computational  advantages  over  QED  with 
regard  to  the  investigation  of  specific  radiative  problems. 
Additionally,  SRED  may  prove  useful  as  a  touchstone  for 
semiclassical  approximations  to  quantum  gravity.16 

In  essence,  SRED  should  be  viewed  as  the  semiclassical 
extension  of  what  has  come  to  be  called  stochastic  electro¬ 
dynamics  (SED).17  In  SED,  one  often  approaches  radia¬ 
tive  problems  through  the  Braffort-Marshall  equation,  a 
generalization  of  the  Abraham -Lorentz  equation  of  clas¬ 
sical  electrodynamics9  that  includes  the  ZPF.  The  result¬ 
ant  theory  thereby  contains  a  fluctuating  vacuum  field 
along  with  the  atom’s  radiation-reaction  field,  and  SED 


has  been  quite  successful  in  accounting  for  various  radia¬ 
tive  processes,  in  particular,  spontaneous  emission  and 
the  Lamb  shift.18  Alternatively,  SRED  is  truly  semiclas¬ 
sical  in  that  the  atom  is  described  quantum  mechanically 
in  terms  of  its  wave  function,  while  the  ZPF  and  the 
radiation-reaction  fields  are  considered  perturbations  on 
the  wave  function’s  evolution.  As  a  consequence,  one  of 
the  main  difficulties  for  SRED  that  differentiates  it  from 
SED  concerns  the  identification  of  a  radiation-reaction  op¬ 
erator  for  inclusion  in  the  Schrodinger  equation. 

In  Section  2  of  this  paper  a  quantum-mechanical  opera¬ 
tor  for  an  atom’s  radiation-reaction  field  is  constructed  by 
use  of  classical  electrodynamics  and  operator  hermiticity 
as  guides.  Then,  in  Section  3,  the  ZPF  is  described  along 
with  its  electric  dipole  interaction  with  an  atomic  system. 
In  Section  4  the  influence  of  the  ZPF  and  the  radiation- 
reaction  field  on  atomic  state  evolution  are  explored  up  to 
second  order  in  perturbation  theory.  One  of  the  major 
problems  with  standard  semiclassical  electrodynamics 
has  always  been  its  inability  to  obtain  spontaneous  decay 
without  spontaneous  excitation  of  the  ground  state.  The 
analysis  of  Section  4  shows  that,  for  an  excited  atomic 
eigenstate  |/n),  the  average  electric  dipole  interactions  as¬ 
sociated  with  the  ZPF  and  the  radiation-reaction  field  add 
to  produce  exponential  decay  of  |m),  whereas  they  cancel 
to  prevent  spontaneous  excitation  of  the  ground  state. 
The  numerical  value  for  the  spontaneous  decay  rate  ob¬ 
tained  in  SRED  equals  that  of  QED,  and  SRED  predicts  a 
radiative  level  shift  for  \m)  in  agreement  with  the  nonrel- 
ativistic  Lamb  shift  of  QED.  The  paper  concludes  with 
Section  5,  in  which  a  numerical  method  for  applying 
SRED  to  more-general  problems  of  radiative  decay  is  dis¬ 
cussed. 


2.  RADIATION  REACTION 

Classically,  when  an  electrical  charge  accelerates  it  emits 
a  field  that  acts  back  on  the  charge.  This  field  is  called 
the  radiation-reaction  field,  ERR,  and  in  nonrelativistic, 
classical  electrodynamics  it  is  simply  related  to  the  rate  at 
which  the  charge’s  acceleration  is  changing: 

2e  2 

Err  =  st3*  =  a<  **-  (1) 

where  fi  is  the  electric  dipole  moment.  As  simple  as  this 
description  seems,  radiation  reaction  presents  many 
thorny  problems  in  classical  electrodynamics.9  For  ex¬ 
ample,  in  the  Abraham-Lorentz  derivation  of  Eq.  (1)  an 
electromagnetic  mass  enters  the  problem  and  becomes  in¬ 
finite  for  a  point  particle  such  as  the  electron.  Addition¬ 
ally,  radiation  reaction  in  classical  electrodynamics  can 
lead  to  preacceleration,  that  is,  acceleration  of  a  charged 
particle  before  the  application  of  any  force.  Here  we 
shall  not  consider  such  problems19  but  rather,  in  the 
spirit  of  Series,7  use  Eq.  (1)  as  a  guide  in  the  construction 
of  an  atomic  operator  whose  expectation  value  is  reminis¬ 
cent  of  the  classical  radiation  reaction  field. 

As  a  first  step  in  transforming  Eq.  ( 1 J  into  such  an  op¬ 
erator,  it  should  be  recognized  that  both  fx  and  need  to 
take  on  atomic  operator  characteristics.  Inasmuch  as  fx 
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is  already  a  quantum-mechanical  operator,  the  task  is  to 
define  an  atomic  operator  to  represent  the  temporal  de¬ 
rivative.  To  this  end,  we  take  advantage  of  the  Schro- 
dinger  equation  for  the  unperturbed  wave  function  and 
the  relationship  that  it  provides  between  <9,  and  the  un¬ 
perturbed  Hamiltonian.  There  is,  however,  some  ambi¬ 
guity  in  this  relationship,  as  we  can  replace  <9,  with  the 
Hamiltonian  as  it  appears  in  the  Schrodinger  equation  for 
the  atomic  wavefunction  or  with  the  Hamiltonian  as  it  ap¬ 
pears  in  the  Schrodinger  equation  for  the  conjugate 
atomic  wave  function: 

iH0 

i ft Bi [ <p)  =  H0\<p)  ^  1  >  yT  >  (2a) 

iH0 

-ifidt{ip\  —  {<p\H0  =>  r)t  •“»  .  (2b) 

A  proper  description  of  spontaneous  emission  in  SRED  re¬ 
quires  the  latter  replacement  [i.e.,  Eq.  (2b)]  because  the 
alternative  choice  leads  to  spontaneous  excitation  of  the 
atomic  ground  state.  Although  it  would  lead  us  too  far 
astray  to  speculate  here  on  the  implications  of  this  choice, 
it  is  intriguing  to  recognize  a  possible  connection  to 
Wheeler’s  and  Feynman’s  absorber  theory,4,20  in  which 
the  seat  of  radiation  reaction  lies  in  fields  that  are  trav¬ 
eling  backward  in  time  from  distant  absorbers. 

Taking  Eq.  (2b)  as  an  operator  equivalent  for  St  and 
maintaining  the  classical  ordering  of  df  and  /*,  one  might 
attempt  to  define  ERR  as  (- 2im3cz)H03ti .  Unfortu¬ 
nately,  Err  is  then  not  Hermitian.  However,  if  this  ex¬ 
pression  for  Err  is  averaged  with  its  Hermitian  conju¬ 
gate,  a  Hermitian  operator  can  be  created.  Thus  we 
define  the  radiation-reaction  field  operator  as 

ERR  =  ~~  ^o3)-  (3) 

For  an  atomic  system  described  by  the  wave  function 
| <p)  the  expectation  value  of  ERR  obtained  from  Eq.  (3)  is 

<<p|Err|?)  =  -rrW  2  [(<p|Ho3|n}(ra|/4|<p) 

Oil  C  n 

-  (<p\v\n)(n\Ho3Ml  (4) 

where  {\n)}  are  eigenstates  of  the  unperturbed  Hamil¬ 
tonian.  If  |<p)  corresponds  to  one  of  those  eigenstates,  for 
example  |m),  then  it  is  straightforward  to  show  that 
(m  |ERR|m)  =  0.  Thus  atoms  in  eigenstates  of  the  unper¬ 
turbed  Hamiltonian  do  not  emit  a  radiation-reaction  field. 
Alternatively,  if  \<p)  corresponds  to  a  superposition  of 
eigenstates  that  can  be  coupled  by  an  electric  dipole  in¬ 
teraction,  for  example  )</?)  —  (l/\-2)[|l)exp(— iwit) 
±  |2)exp(-/w2*)J,  then  under  the  assumption  that 

|lL>2i!  3  6l» 2 W  1  ♦  W21  S  °*2  wi  (5a) 

(e.g.,  |1)  is  the  ground  state  so  that  =  0),  the  expecta¬ 
tion  value  for  the  radiation-reaction  field  becomes 


1  M2ll  ^213 

(<p|ERR|<p)  =  ± — —3 — sin(o>2i^  +  6),  (5b) 

OC 

where  /i2i  =  This  is  essentially  the  classical 

expression  for  an  oscillating  dipole’s  radiation-reaction 
field,  and  it  indicates  that  radiation-reaction  fields  arise 
when  atoms  are  in  superposition  states  of  the  unper¬ 
turbed  Hamiltonian.  Thus  Eq.  (3)  can  be  taken  as  a  rea¬ 
sonable  quantum-mechanical  operator  for  the  radiation- 
reaction  field. 

Although  it  is  to  be  noted  that  Eq.  (5b)  is  half  of  the 
magnitude  of  the  classical  value,  it  should  also  be  recog¬ 
nized  that  Err  is  not  gauge  invariant  because  of  its  de¬ 
pendence  on  H0 .  Thus  the  radiation-reaction  field  op¬ 
erator  of  Eq.  (3)  is  simply  a  calculational  tool,  not  a  true 
physical  quantity.21  Consequently,  any  rigorous  com¬ 
parison  between  the  present  quantum-mechanical  ver¬ 
sion  of  the  radiation-reaction  field  and  its  classical  analog 
must  be  treated  cautiously.22 

For  future  reference  we  note  that  if  VRR  is  defined  as 
the  perturbation  of  the  atom  that  is  due  to  a  dipole  cou¬ 
pling  with  its  radiation-reaction  field;  then  Eq.  (3)  yields 

VRR  =  -fl  ■  Err  =  (n  ■  H03fl  -  M  •  /lH03). 

(6) 

3.  CLASSICAL  ZERO-POINT  FIELD 

As  outlined  by  Boyer,23  Maxwell’s  equations  are  valid  for 
a  variety  of  boundary  conditions.  Standard  classical 
electrodynamics  arises  when  the  boundary  conditions  are 
chosen  such  that  the  fields  in  vacuum  are  zero,  and  only 
retarded  fields  are  allowed.  Random  electrodynamics 
arises  when  the  choice  of  boundary  conditions  consists  of 
random  electromagnetic  radiation  with  a  Lorentz  invari¬ 
ant  spectrum  (i.e.,  the  classical  ZPF),  and  again  only  re¬ 
tarded  fields  are  considered.  Whereas  the  “proper” 
choice  for  boundary  conditions  will  often  lead  to  debate,  it 
should  be  noted  that  according  to  Boyer  “a  consistently 
classical  treatment  of  thermal  radiation  leads  to  the  natu¬ 
ral  introduction  of  temperature-independent  fluctuating 
radiation  in  the  universe.”  24  Here,  questions  regarding 
the  value  of  various  boundary  conditions  for  classical  elec¬ 
trodynamics  will  be  ignored.  We  simply  consider  an 
atom  enclosed  in  a  large  conducting  box  of  length  L,  filled 
with  a  classical  ZPF,  as  a  reasonable  approximation  to  an 
atom  interacting  with  the  vacuum  of  QED.  The  ZPF  vec¬ 
tor  potential  Ao  within  the  box  can  then  be  written  as  a 
modal  expansion  of  random  plane  waves:15 

Ao  =  L'312^  ekx{ckx  exp[i(k  ■  r  -  «,<)] 

kx 

+  ckx*  exp[-i(k  •  r  -  wsO]}.  (7) 

Here  ekx  is  a  modal  polarization  vector,  ckx  is  a  random 
complex  amplitude  for  the  wave,  kt  =  2nt7r/L  {nl 
=  0,  ±1,  ±2, ...),  and 


2 


k  k  = 


2 


[nx 


+  nv 


+ 


n,2] 


(8) 


where  the  modal  index  $  corresponds  to  the  triplet  of 
numbers  {nx  ,  ny  ,  n2}.  Taking  the  derivative  of  Eq.  (7) 
and  replacing  the  sum  over  wave  vectors  by  a  sum  over 
modes,  we  get  for  the  zero-point  electric  field 

1  ^  . 

E0  =  — p=  Z,  exp[i(k„  ■  r  -  «,*)] 

\c2L*  s\ 


or,  equivalently, 

E0  =  i'[(irft)/L3]1/22  esk\f<os[zsK  exp (~iiost) 

sX 

-  zsX*  exp(ia>st)].  (14b) 

Defining  V°  as  the  perturbation  of  the  atom  owing  to  its 
dipole  coupling  with  the  ZPF,  we  have 

V°«)  =  -p  E0  =  -i[(irfi)/i3]1,22  M  ' 

sX 


+  iu)scsk*  exp[-i(ks  •  r  -  wsO]}-  (9) 


The  complex  modal  amplitudes  can  be  written  in  terms 
of  a  scale  factor  csX°  and  a  mean-zero  complex  random 
variable  zsX  :  csX  =  c5X°zsX .  Then,  expressing  zsX  in 
terms  of  two  independent,  unit-variance,  real-valued  ran¬ 
dom  variables  usk  and  vsX  with  mean  zero  such  that 
£{uskuri J  =  t{vskvrft}  =  SsrSXfi,  we  have 

zsx  ~  ( usx  ~~  (10a) 

Z{zSxZrt**}  =  (10b) 

ZbrfrJ  =  =  0,  (IOC) 

where  Z{ }  indicates  an  ensemble  average. 

To  evaluate  the  scale  factor  one  can  use  Eqs.  (9)  and 
(10)  to  show  that 


U  0  = 


1 

4tt 


*{\E0\2}d?r  =  2 

J  s  X 


(<oscsX0)2 

2  7 re2 


(ID 


Lorentz  invariance  of  the  ZPF  spectrum  requires  that  the 
energy  of  each  mode  be  proportional  to  ws.15  Conse¬ 
quently,  the  scale  factor  must  be  inversely  proportional  to 
>/a^,  and  for  agreement  with  the  experimental  results  for 
.blackbody  radiation  the  proportionality  constant  must  be 
related  to  Planck’s  constant  h.  In  this  way  it  can  be 
shown  that 


csX°  =  [(Trhc2)/ws]v2  =>  U0 


(12) 


With  the  aid  of  Eqs.  (10a)  and  (12),  the  electric  field  vector 
of  the  ZPF  now  becomes 


E0  =  [(27 rfi)/L3],/22  cos(ks  '  r  “ 

sX 

-  USK  sin(ks  •  r  -  w,f)].  (13) 

Note  that  this  expression  for  the  ZPF  is  equivalent  to  that 
recently  discussed  by  Ibison  and  Haisch.25 

To  explore  the  interaction  of  the  ZPF  with  a  quantum 
system  we  consider  an  atom  conveniently  placed  at  the 
origin  of  the  coordinate  system  and  make  the  dipole  ap¬ 
proximation  (i.e.,  for  the  modes  of  interest,  k4  *  r  <3  1). 
Then,  in  the  vicinity  of  the  atom, 


E0  =  [(27 rfi)/L3]1/22  cos (<ost) 

S\ 

+  usK  sin(wsi)]  (14a) 


x  \l^s[zsK  exp (-ia>st)  -  zsk*  exp(i«^)]. 

(15) 

Later  we  shall  have  need  for  as  we^  as  f°r 

?{V°mri(t)V°pq(t')},  where  V°mn  is  defined  as  (m|V°|n>. 
Clearly,  from  the  mean-zero  characteristic  of  zsk  , 

%{<*«)}  =  o.  (16) 

Computing  the  second-order  correlation  function  for  the 
perturbation  yields 

=  —73  2  V<*>sa7r(m|/i  '  ««x|n)<Pl/*  •  fr*l<7> 

L  sx 
rn 

X  Z{[z$k  exp {-i(0st)  -  zsX*  exp(itus0] 

X  [zrtl  exp(  —icortr)  -  zrM*  exp(ia>r /')]}. 

(17) 

Multiplying  out  the  terms  in  the  ensemble  average  and 
employing  Eqs.  (10)  yields 

U)Vj,U')} 

=  73-2  \/<<>s a)r{m\/x  •  esk|ra>(p|p-  erJq)- 
L  sx 

rM 

x  {SsrSklt  exp[-i«>s(f  -  <')] 

+  SsrShll  exp[io)s(t  -  C)]}  (18) 

or 

t{V%»  «)Vj,U')} 

irh  _  A  A 

=  77  Z  *>.<»*  I#*  ■  esk|n)(p|p  •  esy\q) 

L  s\ 

x  {exp[-ik>s(t  -  t ')]  +  exp [iois{t  -  *')]}•  (19) 

4.  PERTURBATION  THEORY 

A.  Evolution  of  the  Wave  Function 

We  now  consider  the  Schrodinger  equation  for  an  atom  in 
the  presence  of  the  ZPF  and  the  radiation-reaction  field, 

ifi^-  =  [Ho+  v™+  v°(/)]k>,  (20) 

dt 

and  expand  the  atom’s  wave  function  in  terms  of  the  un¬ 
perturbed  basis  wave  functions: 
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W)  -  2  a«(Ol«)exp(“/aiB0.  (21) 

n 

Then  in  standard  fashion  we  obtain  an  equation  for  the 
expansion  coefficients: 

a"(t)  =  --  2  [(n|VRR|,/)  +  <n|V0|jf)]a/(f)  exp (iwnJt), 

(22) 

where  wnj  =  u>n  -  <oj .  Integrating  and  defining  an°  as 
an(0),  we  can  also  write  Eq.  (22)  as 

o.«)  =  on°  -  lr  2  f'[v™  +  vj/nwn 

n  j  Jo 

X  exp(iu)njt’)dt' .  (23) 

B.  Wave-Function  Evolution  Owing  to  Radiation 
Reaction 

Defining  the  first  term  on  the  right-hand  side  of  Eq.  (22) 
as  anRR(f ),  we  have,  after  substituting  from  Eq.  (6), 

X  (rc|(/t  •  tf0V  -  A«  •  t*Ho3)\j)  exp (ia>njt). 

(24) 

Note  that  the  radiation-reaction  field  forces  an  orienta¬ 
tion  onto  the  atom’s  dipole  moment.  Taking  an  ensemble 
average  of  Eq.  (24)  over  the  ZPF,  and  defining  bn(t)  as 
^{an(0},  yields 

i.“(0  =  rrrr  2  bj(t)[finp(p\H03\q)nqJ 

ofl  C  pq  j 

~  MnpMp,(9l^o3li>]exp(ia>„/)  (25) 

or,  when  the  Hamiltonian  terms  are  evaluated, 

*„RR(0  =  ^5  2  bj(t)finp/j.pj(wp3  -  a >f)  exp  (ia>nJt). 

(26) 

We  now  ignore  the  rapidly  rotating  terms  in  Eq.  (26)  [i.e., 
exptia^  O  — +  SnJ  ]  and  recognize  that  the  sum  in  Eq.  (26) 
will  be  dominated  by  those  states  |jd)  that  best  satisfy  the 
validity  condition  for  the  radiation-reaction  field  operator, 
inequality  (5a).  Thus  we  approximate  ( o)p 3  -  a^3)  with 
< opn 3  and  obtain 

*nRR(0  =  -^3  2  U/ip|2"np3’  (27) 

C.  Wave-Function  Evolution  Owing  To  the  Zero-Point 
Field 

We  now  define  the  second  term  on  the  right-hand  side  of 
Eq.  (22)  as  anz??{tY. 

anZPF  ~  --  2  V°nJ(t)aj(t)  exp(iunJt),  (28) 

n  j 


and  substituting  from  Eq.  (23)  yields 


amZPF(t)  =  --  2  K,n(t)  exp (i<omnt) 

*  n 


a„°-  ^2  fV“  +  Vt/t'W) 

h  j  Jo 


J  Jo 
X  exp(i<ony-f ')df ' 


(29) 


We  can  now  ensemble  average  Eq.  (29)  over  the  ZPF,  so 
that 

i>mZPF(t)  =  -Tj  2  exp(ia.m„() 

*  nj 

X  fV{V®,(D V°nj(t')}bj(t')  exp (iwnjt')dt'. 

Jo 

(30) 

In  writing  Eq.  (30)  we  have  taken  advantage  of  the  fact 
that  causality  requires  a/t')  to  be  independent  ofV^Jt) 
for  t '  <  t ,  and  we  have  assumed  that  the  evolution  of  the 
ZPF  is  unrelated  to  the  atom’s  dynamics.  Thus  we  con¬ 
clude  that  a  jit')  and  V°mn(t)  are  independent  of  each 
other  for  V  <  t.  Moreover,  inasmuch  as  the  ZPF  pertur¬ 
bation  is  weak  in  the  limit  of  a  large  volume  enclosure,  to 
second  order  in  perturbation  theory  we  can  make  a  deco¬ 
rrelation  approximation  between  V®j(t')  and  ajit’)  {i.e., 
the  ZPF  correlation  time  is  much  less  than  h[V°nj{t  ')]-1}. 
In  this  way  we  obtain 

X[V°mn(t)V0nj(f)aj(t')}  ^  t[Vl 

=  ■avln(t)Vanj(f)}bj(t').  (31) 

To  second  order  in  perturbation  theory  we  approximate 
bj(t ')  with  its  value  at  t  and  substitute  for  the  correlation 
function  from  Eq.  (19).  Equation  (30)  then  becomes 

&mZPF(f)  =  —  ““3  2  bj{t)  exp (iojmnt) 

X  I  (2  '  esX|n>(n|A»  •  ej» 

X  {exp[-iws(*  -  tf)]  +  expfi&^U  -  t')]}j 

x  exp(i(Dnjt')dt'.  (32) 

As  we  discussed  above,  the  dipole  moment  of  the  atom 
is  aligned  along  the  direction  of  the  radiation-reaction 
field.  We  define  this  orientation  as  the  z  axis  of  the  co¬ 
ordinate  system26: 

^  k* 

2  =  1  ”  j^|2»  (33) 

so  Eq.  (32)  becomes 
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bmz?Y(t )  =  — —3  X  bA 0  exp(io>m„i)MmnMn> 

flL  nj 

x  J  |X  w»  sin2(0s){exp[-iws(<  -  i')] 

+  exp[io>,(t  -  <')]}j  exp(ia»n/')dt',  (34) 

where  6S  is  the  polar  angle  between  the  wave  vector  for 
the  sth  mode  and  the  z  axis  as  defined  by  the  radiation- 
reaction  field.  If  we  now  allow  the  volume  of  the  box  en- 
closing  the  atom  to  go  to  infinity, 

Eq.  (34)  becomes 


to3  exp( -iu)t) 


Int(l)  =  lim  [ 
a— 0  J  0 

X  [  'exp[a/'  +  i{(x)  +  wnJ)r]drdw,  (39) 

Jo 

which  on  integration  with  respect  to  the  temporal  vari¬ 
able  yields 

I  fx  o>3doj 

Int(l)  =  limlexp{[(a  +  — 


+  i(w  +  (onj) 


r. 


CO3  exp(  -i(0t)da> 


(40) 


J0  a  +  i(a)  +  (*)nj )  J 

Note  that  both  integrals  will  be  dominated  by  large  w. 
However,  in  the  second  integral  these  terms  oscillate  rap¬ 
idly  and  quickly  average  to  zero.  Thus 


i>mZPF(t)  =  ~g7r2^3  2.  bj(t)  exp (i(OmJ)^mn^nj  Int(  1 )  =  lim|exp{[(a  +  iwnj)t]} 


-o\ 


X  J  y  J  J  w3{exp[-ifc>(*  -  *')] 
+  exp[iw(£  —  £')]}  sin3(0)d0d</>d< 


r*  w3[a  -  i(u)  +  wnj)]do> 

Jo  a2  +  (w  +  <onj)2 


(41) 


X  ex.Tp(iconjt’)dt' . 

Integrating  over  8  and  <j>  then  yields 
1 


(36) 


6mZPF(<)  =  -——3  X  &,■(*)  exp (i(omnt)fimnfi-nj 
3  77/1 C  nj 

X  J  (J  w3{exp[-ta>(*  “  *')] 

+  expfiwU  -  *')]}da>)  exp(i(onjt')dt’ 


As  discussed  by  Heitler,27  the  limit  of  the  real  term  yields 
a  Dirac  delta  function,  whereas  that  of  the  imaginary 
term  can  be  expressed  in  terms  of  a  principal  value: 

Int(l)  =  exp(iwn/)|7r J  oq3S(u)  +  wnj,)da> 

a>3dw 


-‘4"« 


O)  +  (*in  j) 


(42) 


(37) 


or,  when  we  rearrange  the  orders  of  integration, 

1  ^ 

bmZFF(t)  =  --  -—3  bj(t)  exv(i(t)mnt)fimnHnj 
3  777?  C  nj 

X  J  a>3jexp(  -io)t)  J  exp[f(w  +  wn;*)f']dJ' 
+  exp(ia>oJ  exp[-i(w  -  wn,-)*']d*' j da>, 


In  a  similar  fashion  the  second  integral  is  evaluated  and 
yields 

Int(2)  -  exp(iwn^o|^ J  u)3  8(a)  -  wn^)do> 
i>3dfo 


+  i'jP 


^  f°°  id3 dio 

iy,)Q  (Z  -  z~) 


(43) 


Combining  Eqs.  (42)  and  (43)  in  Eq.  (38b),  we  get 
1 


bmZPT(t) 


„■  X  bj(t)  exp [iiomjt)ii  mn  Vnj 
She  nj 

J  cu3[<5(w  +  (onj )  +  8(a)  -  o>n7)]df/> 

Jo 


(38a) 


bmZPF(t)  =  -T-T-3 


X  bj(t) 


3  77 he 3  nj 

X  exip(ioi)mrit)fimnfinj[lnt(l)  +  Int(2)]. 

(38b) 

To  proceed,  we  now  evaluate  the  integrals  in  Eqs.  (38) 
separately.  Considering  the  first  integral,  we  write  this 
as 


-  i A, 

where  the  imaginary  term  is  defined  as 
1 


(44) 


A  = 


Sirhc3 


2  bj(t)  exp (i  Mm  jt)  fimnV'nj 


X  jA 


i; 


*>nj  w  w nj 


dw 


(45) 
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and  leads  to  a  radiative  level  shift  induced  by  the  ZPF. 

Recognizing  that  the  rapidly  oscillating  terms  in  Eqs. 
(44)  and  (45)  quickly  average  to  zero,  we  have 


bmZPF(t)  =  - 


bm(t) 


s  u..i  2[  r 

n  [JO 


u)3S(od  +  <o„m)dw 


3  he3 

+  |  cj3S(u>  -  Wfjmldw  —  (46a) 

Jo 


where 


A'  e 


37 rhca  n 


2  f 

n  JO 


w3do> 

2  2  ' 
CO 


(46b) 


Note  that  in  Eq.  (46a)  the  integrals  are  only  over  positive 
values  of  <d.  Thus,  depending  on  the  energetic  relation¬ 
ship  between  |n)  and  |m),  the  delta  functions  will  cause 
only  one  of  the  two  integrals  to  be  nonzero.  Thus 


imZPF«)  = 


MO 

Me3  r 


I  Mi* 


3 

(t)  3 

wnm 


com  >  OJn 

a*n  > 


-  ibm(t) A\ 


(47) 


D.  Spontaneous  Decay 

Combining  Eqs.  (47)  and  (27),  and  ignoring  the  radiative 
level  shift  term  for  the  moment  as  this  has  no  influence  on 
decay,  we  find  for  the  evolution  of  the  eigenfunction  coef¬ 
ficients  that 


*«<o  =  - 


bjt) 

She3 


X  2  lMmni|  ^ 


<*>n  >  (x)n 


(48a) 


2  bm(t) 

bm(  0  —  3  I  Pmn 

one  n<m 


2 . ,  3 


(48b) 


Defining  Pm  as  the  probability  for  finding  the  atom  in 
the  excited  state  | m)  (i.e.,  |6m|2),  we  obtain  an  exponential 
decay  equation  for  this  probability: 

=  -r  Pm  ,  (49a) 

d<  m  m  » 


where 


Tm  -  2  \PmnW-  (49) 

This,  of  course,  is  just  the  Einstein  A  coefficient  for  the 
spontaneous  decay  rate  of  |m).  Note  that  if  we  define  |1) 
as  the  atomic  ground  state,  then  T  j  is  identically  zero  be¬ 
cause  there  are  no  states  of  lower  energy.  Essentially,  in 
SRED  excited  states  decay  exponentially  because  the  ra¬ 
diation  reaction  and  vacuum  pathways  add  construc¬ 


tively.  Spontaneous  excitation  of  the  ground  state  does 
not  occur,  because  in  that  case  these  two  pathways  would 
interfere  destructively.  Although  many  authors  have 
recognized  this  interference  between  radiation  reaction 
and  the  vacuum  in  QED,  here  we  find  it  again  in  the  con¬ 
text  of  a  semiclassical  description  of  spontaneous  decay. 

E.  Lamb  Shift 

Returning  now  to  the  radiative  level  shift  described  by 
Eq.  (46b),  we  note  that  the  integral  is  strongly  influenced 
by  the  very  high-frequency  ZPF  modes: 

(50) 

37 rhe6  n  Jo 

Employing  the  Thomas-Reiche-Kuhn  sum  rule28  and  re¬ 
membering  that  \fjLmn\  —  e\zmn\y  we  can  rewrite  relation 


This,  however,  is  just  the  average  value  of  the  interaction 
energy  between  a  free  electron  and  the  z  component  of  the 
ZPF,  as  can  easily  be  shown  from  Eqs.  (7)  and  (33)  and 
expression  (35): 


“•-(^)<ia”,v  <5a 

Because  this  energy  is  the  same  for  all  levels,  it  does 
not  give  rise  to  an  observable  shift  and  so  should  be  sub¬ 
tracted  from  the  shift  of  Eq.  (46b).  Of  course  this  sub¬ 
traction  is  just  the  well-known  mass  renormalization  of 
QED.12’29 

The  observable  radiative  level  shift  in  SRED  is  thus 


A  - 


3  irhcd 


3K 


„l  V  f 

Jo 


codco 


,  (53) 


and,  when  the  integral  is  evaluated  with  a  high-frequency 
cutoff  of  mc2/h,  the  result  is  the  standard  Lamb-shift  for¬ 
mula  of  nonrelativistic  QED: 


A  - 


2 

37rftc3 


2 


nl2  M 


me2  \ 

\En~^E~\j- 


(54) 


In  SRED,  as  can  be  argued  in  QED,12  the  Lamb  shift  is 
seen  as  an  ac  Stark  shift  of  an  atomic-energy  level  in¬ 
duced  by  the  ZPF.30 


5.  NUMERICAL  SIMULATION  OF 
SPONTANEOUS  DECAY 

One  of  the  attractive  features  of  SRED  with  regard  to 
computation  concerns  the  representation  of  the  vacuum 
as  an  easily  simulated  classical  ZPF.  A  potential  utility 
of  SRED  may  therefore  reside  in  the  numerical  simula¬ 
tion  of  spontaneous  decay  in  an  arbitrary  vacuum  (e.g.,  a 
cavity).31  Unfortunately,  the  description  of  spontaneous 
decay  in  SRED  as  formulated  above  requires  an  expres- 
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sion  for  radiation  reaction,  and  in  an  arbitrary  vacuum  it 
may  be  difficult  to  arrive  at  a  workable  operator  represen¬ 
tation  of  radiation  reaction  for  inclusion  in  the  Schro- 
dinger  equation.  To  circumvent  this  problem,  in  this  sec¬ 
tion  we  appeal  to  the  fluctuation-dissipation  theorem  in 
order  to  remove  the  explicit  reference  to  radiation- 
reaction  from  the  equations  that  describe  a  wave  func¬ 
tion’s  evolution.  As  will  be  shown  below,  for  the  vacuum 
of  free  space  the  procedure  provides  an  accurate  numeri¬ 
cal  simulation  of  the  sodium  S2Pm  state’s  radiative  de¬ 
cay. 

For  simplicity  we  consider  the  case  of  a  two-level  atom, 
where  El  =  0  and  E2-  ha)  21  >  and  we  take  advantage  of 
the  wave  function’s  normalization  to  write  |cr  1 12  —  1 
-  |a2|2.  Thus  with  regard  to  spontaneous  decay  we 
need  only  consider  a  second-order  perturbation  solution  to 


a2  =  a2RR  +  a2ZPF.  (55) 

As  the  previous  analysis  demonstrates,  if  Eq.  (55)  were  to 
be  averaged  over  the  vacuum  field  then  each  term  on  the 
right-hand-side  would  contribute  equally  to  the  decay  of 
a 2;  that  is, 

Re(a2KR)  =  Re(o2ZPF)  =  ~  <o2«)>-  (56) 


As  many  authors  have  stated,  this  equality  is  to  be  inter¬ 
preted  as  a  manifestation  of  the  fluctuation-dissipation 
theorem;32  that  is,  the  wave  function’s  evolution  as  a  re¬ 
sult  of  radiation  reaction  corresponds  to  the  atom’s  dissi¬ 
pative  response  to  the  random  fluctuating  force  of  the 
vacuum  field.  Physically,  then,  Eq.  (56)  would  suggest 
that  radiation  reaction  can  be  eliminated  from  perturba¬ 
tion  theory  by  the  replacement  of  terms  that  correspond 
to  radiation  reaction  with  terms  that  correspond  to  the 
vacuum  field,  so  long  as  there  is  a  proper  accounting  of 
each  term’s  relative  sign.  Following  this  procedure  for 
Eq.  (55),  we  obtain  a2  =  2a2ZPF  or,  using  Eq.  (29), 

2  i 

a 2(f)  =  y  m216*  *  Eo(f)  exp(za>21*) 


x 


+ 


[Ml2£z  '  E0(f')  V'Jf] 


X  a2(t')  exp(-i(i)2it,)dt' 


(57) 


As  our  interest  is  in  the  atom’s  spontaneous  decay,  we 
can  perform  a  coarse-grain  averaging  of  Eq.  (57)  over  a 
time  scale  r  that  is  long  compared  with  the  correlation 
time  of  the  vacuum  field  but  short  compared  with  the 
wave  function’s  evolution.  Then,  with 

l  ft  +  r/2 

{Y{t))r  =  -  Y(t')dt',  (58) 
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Eq.  (57)  becomes 

2[yU.2i|2  ft  *  * 

a2{t)  = - 7n  {exp{i(o21t)ez  •  E 0(t)et  ■  E0(*') 

n  Jo 

X  exp{-i(D2\t'))Ta2(t ')d£\  (59) 


Substituting  for  ez  •  E0U)  from  Eq.  (15),  defining  Ss  as 
<*>s  -  o)2l ,  and  ignoring  the  rapidly  rotating  terms  of  ws 
+  o>2i ,  we  obtain 


a2(t)  =  - 


2ir|/i2il2  v 
*L3  sKpv 


Zs\Zpv*  Vwsa>p(esX.2ep,,_.) 


.  (60) 


Because  of  the  independent  nature  of  the  zsX  ,  as  the 
vacuum  field  mode  spacing  Aco  decreases  we  expect  that 


lim  zsX2  *P,*  =  (61) 

Aw->0  P  V 


Thus  the  sum  over  vacuum  field  modes  in  Eq.  (60)  can  be 
approximated  with 


a2{t) 


2tt|m2iI2 

hL3 


2 


X 


Jo 


(exp [-iS8{t  -  t’)'\)Ta2(t,)dt' 


.  (62) 


Making  the  second-order  perturbation  theory  replace¬ 
ment  of  a2{t')  with  a2(t)  on  the  right-hand  side  of  rela¬ 
tion  (62)  and  evaluating  the  integral  then  yield 


a2(t) 


27T\fl21\2 

f\L3 


a2(t >2 

sX 


(6SX,2€SX, 


.) 


x  [1  -  exp(-iSst)] 


(63) 


For  the  purposes  of  computational  convenience  we  now 
set  |2:si|2  “  |zs2|2,  recognizing  that  doing  so  should  have 
little  effect  on  a  numerical  simulation  of  spontaneous  de¬ 
cay  so  long  as  the  number  of  modes  used  in  the  computa¬ 
tion  is  large.  The  sum  of  polarization  terms  in  Eq.  (63) 
can  now  be  easily  handled  as  in  Eq.  (33)  and  then  aver¬ 
aged  over  4tt  steradians,  so  that 


a2{t)  = 


4w|^2il2 
3  hL3 


a2(*)X 
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[1  -  exp(-i’40] 


(64) 


We  can  further  simplify  Eq.  (64)  by  recognizing  that  the 
volume  of  the  vacuum  under  consideration,  L3,  is  propor¬ 
tional  to  the  vacuum  field  mode  spacing  (i.e.,  L3 
=  27t2c3/<d2Acix),33  so  that 

2\ii21\2m212\<x)  I 
«2(0  = - ^^3 - CL2{t)\\z0\Zm2lt 


+  2' 


[1  -  exp(-z«5s*)] 


(65) 


where  the  prime  on  the  summation  indicates  that  the  Ss 
=  0  term  has  been  removed. 

Finally,  integrating  Eq.  (65)  yields 
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In 


0=2(0 

a2(0). 

—  K 

l*ol! 


d>2\t 


+  2'  w2 W^P 


where 


-  cos (Sst)]  +  — 2  [sin(V)  ~  4*] 


2|M2l|2w212Aa> 

37T^C3 


(66) 


(67) 


Again,  with  P2(t)  defined  as  the  probability  of  finding  the 
atom  in  the  excited  state,  Eq.  (66)  results  in 

P2(t)  =  P2(0)exp|-K|2ol2^2i*2  ~  2*2  *  |l**l2(j2 


X  [1  -  cos(£«0] 


(68) 


To  test  the  validity  of  this  semiclassical  procedure  for 
numerically  simulating  spontaneous  emission  we  can  ap¬ 
ply  Eq.  (68)  to  the  decay  of  sodium’s  first  excited  state  at 
589.6  nm  (i.e.,  3 2Pm  —  3  2S„2).  The  3 2Pin  state  life¬ 
time  is  17.0  ns,34  yielding  an  electric  dipole  moment  for 
this  transition  of  2.07  X  10-29  Cm.35  Figure  1  shows 
P2{t)  computed  from  Eq.  (68)  for  three  values  of  the  mode 
spacing  (i.e,  100,  10,  and  1  kHz).36  The  range  of  modes 
employed  in  the  summation  of  Eq.  (68)  was  limited  to  ±10 
GHz,  which  corresponds  to  a  correlation  time  for  the 
simulated  vacuum  field  of  —0.01  ns.  Clearly,  as  the 
mode  spacing  approaches  zero  the  decay  ofP2(0  becomes 
exponential  to  high  precision,  and  the  numerically  simu¬ 
lated  rate  of  decay  corresponds  to  a  17-ns  lifetime. 


6.  SUMMARY 

The  Einstein  A  coefficient  for  spontaneous  decay  has  been 
derived  by  use  of  second-order  perturbation  theory  in  a 
completely  semiclassical  formalism  with  (1)  the  assump¬ 
tion  of  a  classical  ZPF  and  (2)  a  Hermitian  operator  for 
the  atom’s  radiation-reaction  field.  Additionally,  the 
semiclassical  theory  has  yielded  the  Lamb  shift,  which 
appears  as  an  ac  Stark  shift  induced  by  the  ZPF.  Al¬ 
though  these  results  are  intriguing  in  their  own  right,  the 
primary  motivation  for  developing  SRED  is  for  its  use  as 
a  tool  in  computer  simulations  of  spontaneous  decay  un¬ 
der  various  conditions.  It  was  therefore  demonstrated 
that  SRED  has  the  potential  to  simulate  spontaneous  de¬ 
cay  in  arbitrary  vacuum  fields  once  the  modal  composi¬ 
tion  of  the  vacuum  field  is  specified. 
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TECHNOLOGY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer”  for  national  security  programs,  spe¬ 
cializing  in  advanced  military  space  systems.  The  Corporation's  Technology  Operations  supports  the 
effective  and  timely  development  and  operation  of  national  security  systems  through  scientific  research 
and  the  application  of  advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical 
staffs  wide-ranging  expertise  and  its  ability  to  stay  abreast  of  new  technological  developments  and 
program  support  issues  associated  with  rapidly  evolving  space  systems.  Contributing  capabilities  are 
provided  by  these  individual  Technology  Centers: 

Electronics  Technology  Center:  Microelectronics,  VLSI  reliability,  failure  analysis, 
solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared  and 
CCD  detector  devices,  Micro-Electro-Mechanical  Systems  (MEMS),  and  data  storage 
and  display  technologies;  lasers  and  electro-optics,  solid  state  laser  design,  micro-optics, 
optical  communications,  and  fiber  optic  sensors;  atomic  frequency  standards,  applied 
laser  spectroscopy,  laser  chemistry,  atmospheric  propagation  and  beam  control, 
LIDAR/LADAR  remote  sensing;  solar  cell  and  array  testing  and  evaluation,  battery 
electrochemistry,  battery  testing  and  evaluation. 

Mechanics  and  Materials  Technology  Center:  Evaluation  and  characterization  of  new 
materials:  metals,  alloys,  ceramics,  polymers  and  composites;  development  and  analysis 
of  advanced  materials  processing  and  deposition  techniques;  nondestructive  evaluation, 
component  failure  analysis  and  reliability;  fracture  mechanics  and  stress  corrosion;  analy¬ 
sis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch  vehicle 
fluid  mechanics,  heat  transfer  and  flight  dynamics;  aerothermodynamics;  chemical  and 
electric  propulsion;  environmental  chemistry;  combustion  processes;  spacecraft  structural 
mechanics,  space  environment  effects  on  materials,  hardening  and  vulnerability  assess¬ 
ment;  contamination,  thermal  and  structural  control;  lubrication  and  surface  phenomena; 
microengineering  technology  and  microinstrument  development. 

Space  and  Environment  Technology  Center:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and 
ionospheric  physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing, 
hyperspectral  imagery;  solar  physics,  infrared  astronomy,  infrared  signature  analysis; 
effects  of  solar  activity,  magnetic  storms  and  nuclear  explosions  on  the  earth’s  atmos¬ 
phere,  ionosphere  and  magnetosphere;  effects  of  electromagnetic  and  particulate  radia¬ 
tions  on  space  systems;  component  testing,  space  instrumentation;  environmental  moni¬ 
toring,  trace  detection;  atmospheric  chemical  reactions,  atmospheric  optics,  light  scatter¬ 
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